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Abstract 

Because of the rotational invariance of the nuclear Hamiltonian, there 
exists a density functional (DF) for nuclei that depends only on two scalar 
densities. Practical calculations boil down to radial, one-dimensional ones. 

Let Z, N, A = Z + N be the proton, neutron and mass numbers, respectively. 
The nuclear Hamiltonian H is known to be invariant under rotations. Therefore, 
besides Z and N, nuclear ground states (GSs) carry good quantum numbers, 
J and M, for the total angular momentum and its z-component. Two cases 
occur : either J = 0, hence the GS is not degenerate and its density is isotropic, 
or J > 0, hence one sees a trivial degeneracy for a magnetic multiplet of GSs, 
the densities of which, non isotropic, contain multipole components pQ [5] up to 
order 2 J, with the same monopole for all members of the multiplet. 

In any case, the ensemble density operator for the GS(s) of a given nucleus, 

B ZN = (2J +1)- 1 ^2\ZNJM)(ZNJM\, (1) 

M 

is a scalar under rotations. In any case also, the GS energy obtains by, 

E ZN = Tr B ZN H. (2) 

Every spherical harmonic function, except a monopole, integrates out to 
zero. Therefore, only the monopole components of the proton and neutron 
densities contribute to the density normalizations which identify a nucleus, 

/>OC 

r 2 dr p pQ (r) = Z, / r 2 dr p n0 (r) — N, (3) 
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the position r being here taken as a scalar rather than a vector. In the following, 
for conciseness, we omit the subscript and often denote p the pair {p p , p n }. 
Consider the "density constrained search" [3] [4] for a minimal energy, 



Inf p ^{z,N} 
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( Infg^ p Tr B H ) + / r 2 dr [u p (r) p p (r) + u n (r) p n (r)' 
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(4) 
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where it is understood that the many-body density operator B is restricted to 
be a scalar under rotations. Also u p and u n are scalar. With two distinct 
densities p v ,p n , two potentials are needed to represent the external potential 
used by the Hohenberg-Kohn [5] theorem as a functional Lagrange multiplier 
to constrain the density. Notice that now we do not put subscripts Z, N to 
B, because the normalizations, Eqs. ([3]), are implemented at the stage of the 
"outer" minimization. At the "inner" stage, TV and Z do not need to be integers. 
This inner minimization, Infg^ p , defines a density functional, 

F[p] ee Inf H ^ p Tr B H, (5) 

and, in that sector defined by the additional constraints, Eqs. (j3|), with now Z 
and N physical integers, the GS energy of a nucleus results from, 

/>oo />oc 

E ZN =Irf p F[p], / r 2 drp p (r) = Z, / r 2 dr p n (r) = N. (6) 
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It seems, therefore, that Eqs. §5§ and ([6]) provide the basis of a fulfledged nuclear 
DF theory (DFT) in radial space, a symmetrized theory d la Gorling [B]. 

However, one must first remove an ambiguity in the definition of the degree 
of freedom r and the associated density p(r). For A nucleons, the simplest set 
of degrees of freedom are the single nucleon coordinates fi, fa, ... ta and the 
simplest definition of the density consists in integrating out all of them but one, 

p(r) — A J dridr 2 ...drA-i |*(fi, fa, fk-i, r)| 2 ■ (7) 



But, as discussed by [7], it is more physical to use a density a \f — Rj , measured 

from the center-of-mass (CM) coordinate R = (fi + ... + ta)/A of the nucleus, 
rather than the density p(r), defined in the laboratory frame. Since H is also 
translation invariant, the wave function for A nucleons is rather an "internal" 
one, ipint •••) £4.-1)5 °f iA ~ 1) Jacobi coordinates only, £1 = r% — fi, £2 = 
r 3 - (fa + fi)/2, ... , U-i =r A - (f A -i + ?a-2 + ••■ + n)/{A - 1) . ^ 

It turns out that the last Jacobi coordinate is proportional to ?a — R, namely, 
^A-i = a^T ~ R)- Except for trivial scaling factors, the "internal" density 
appears naturally to be, 

a(i) = Aj d£l dg 2 ... d£ A -2 \Ant(£l, 6, 6-2, Ol 2 ■ (8) 

Strictly speaking, what we called F[p] should rather be a functional of a. Un- 
fortunately, there is no need to stress that calculations with Jacobi coordinates 
are much more complicated than calculations in the laboratory frame. 

The solution found in 8 to link internal density a and laboratory density p 
consists in trapping the CM by a harmonic potential. The Hamiltonian becomes, 

»=1 \i=l / i>j=l i>j>k=l 
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Here v is the usual two-body interaction, and one can also include a three- 
body interaction w or even the luxury of more-body ones. Except for the trap, 
Galilean invariance is requested; no density dependence in v,w, ... is allowed. 
The trap is parametrized by an arbitrary, but fixed constant K, to be chosen 

/ A \ 2 

for maximum convenience of practical calculations. The term, Km I Xh=i ) > 

can also be written as, Amui 2 R 2 , with the CM coordinate R = A" 1 Ylt=i ft an d 
lu = \J K A; the CM frequency depends on the mass number A. The same form 
also shows that the trap is the sum of a one-body and a two-body operators, 
representable in second quantization without coefficients depending on A, 



K m 



E 5 



Km 



A 

E< 

i>j=i 



(10) 



In each sector specified by integer Z and N, the GS(s) of Ti factorize(s) as 
product(s) of a common Gaussian T for the CM and internal wave function(s) 
of the (A — 1) Jacobi coordinates, 



V(r 1 ,...,f A )=T(R) ^i„ t (&,...,6t-i), T(R) =tt-3&-3 ex p 



2b 2 



(11) 



with b — [h/(A m to) ] 5 . The Gaussian is rotation invariant. Its does not perturb 
the physical quantum numbers J, M when H. is substituted for H. 

As shown in [5j the link between p and a is a trivial, invertible convolution, 



p(f) 



,4 3 



(.4-1) 



<iR [r(i?)] 2 cr 
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A- 1 



i? 



(12) 



This link is the same for any member of a magnetic multiplet. Hence it extends 
to the scalar densities p and a provided by that kind of ensemble described by 
Eq. {1} . The vector coordinates we used temporarily to handle the non isotropic 
densities of individual members of a magnetic multiplet can be reduced again 
to scalars. The convolution, p = T 2 * a, see Eq. (fT2|) . actually becomes, 



■p(r) 



27T-5A 3 
(A - l) 3 b 



ds exp 



r 2 + s 2 



sinh 



2rs 



A 



A-l 



(13) 



Since it is easy to second quantize Ti. and use the laboratory proton and neu- 
tron densities, p p (r) — Tr B a pr c pr , p n (r) = Ti B cj^ c nr , where, with obvious 
notations, we have introduced proton and neutron creation and annihilation 
operators at a scalar position r, traces upon scalar density operators can be 
calculated in Fock space. Then one only tunes Eqs. ([6]) into, 



E ZN + ^hVKA = In£ p J r [p], / r 2 drp p {r) = Z, / r 2 dr 
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p n (r)=N. (14) 



We conclude by claiming that a DFT for nuclei energies (and every other 
scalar) is available with scalar densities. Practical designs of a nuclear DF, 
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F[p] = Infg^pTr BH, simultaneously valid for doubly even, odd and doubly 
odd nuclei, can be attempted in the laboratory radial frame. 

From a DFT, one expects densities and energies. At the cost of ignoring 
density multipole^l other than monopoles, our approach does provide formally 
exact energies, in a most simplified, one-dimensional theory. A major problem 
of nuclear physics is the prediction of exotic nuclei. Our result is of a special 
interest for the study of the neutron drip line, where neutron halos are notori- 
ously difficult to describe. Dimensional reduction should make it easier to better 
focus the theory on the desigrjl of the functional J 7 , a still formidable problem. 
Whether the same reduction to "radial pictures" might also simplify the non 
local [TO] versions of the DFT, in particular the quasi-local versions derived from 
Skyrme force models and labelled "energy density" theories, is likely 
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lr The wave function multipoles, nonetheless, contribute to the density monopole. 

2 See [9] for expansions of densities in polynomials constrained by matter conservation. 
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